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Abstract 

We predict the amount of cometary, interplanetary, and interstellar cosmic dust 
that is to be measured by the Cometary and Interstellar Dust Analyzer (CIDA) 
and the aerogel collector on-board the Stardust spacecraft during its fly-by of 
comet P/Wild 2 and during the interplanetary cruise phase. We give the dust 
flux on the spacecraft during the encounter with the comet using both, a radi- 
ally symmetric and an axially symmetric coma model. At closest approach, we 
predict a total dust flux of 10^ ° m^^ s^^ for the radially symmetric case and 
2^06.5 ^-2 g-i axially symmetric case. This prediction is based on an 

observation of the comet at a heliocentric distance of 1.7 AU. We reproduce 
the measurements of the Giotto and VEGA missions to comet P /Halley using 
the same model as for the Stardust predictions. The planned measurements of 
interstellar dust by Stardust have been triggered by the discovery of interstellar 
dust impacts in the data collected by the Ulysses and Galileo dust detector. 
Using the Ulysses and Galileo measurements we predict that 25 interstellar par- 
ticles, mainly with masses of about 10"^^ g, wiU hit the target of the CIDA 
experiment. The interstellar side of the aerogel collector will contain 120 in- 
terstellar particles, 40 of which with sizes greater than 1 fim. We furthermore 
investigate the "contamination" of the CIDA and collector measurements by 
interplanetary particles during the cruise phase. 



1 Introduction 

The Stardust mission to comet 81P/Wild 2 is dedicated to the in situ measure- 
ment and sample return of cosmic dust. At the comet the pristine cometary 
material is investigated which is assumed to be a good sample of the nebula 
from which the Solar System has formed. Pre-solar interstellar grains that sur- 
vived the formation process may therefore be present in the samples collected 
or measured in situ. On the way to the comet the opportunity is taken to 
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investigate contemporary interstellar dust that traverses the Solar System and 
originates from the local interstellar low-density diffuse cloud which crosses the 
Sun's way through the Milky Way We predict the dust flux on the two 

main instruments on-board Stardust, the aerogel dust collector and the Cos- 
mic and Interstellar Dust Analyzer (CIDA) [Q. The dust collector is a plate 
with an area of « 0.1 m^ that is to be exposed to the dust stream. Particles 
hitting the collector are trapped inside the aerogel CIDA is an in-situ im- 
pact plasma dust detector with a time-of-flight mass spectrometer similar to the 
PUMA/PI A instruments flown on Giotto Stardust launches in February 

1999 and flics by Earth in January 2001 for gravity assist. The slow (6.1 km s~^) 
fly-by of comet P/Wild 2 occurs on January 1st, 2004 at a heliocentric distance 
of 1.86 AU when the comet is on the out-bound part of its orbit. During the 
fly-by the distance of the spacecraft to the comet nucleus is ~ 100 km, and the 
phase angle, i.e. the Sun-comet-spacecraft angle, is 70°. After fly- by Stardust 
returns to Earth in January 2006, where the sample-return capsule that con- 
tains the aerogel collector is separated from the spacecraft and sent to a direct 
entry into Earth's atmosphere. The collection and in situ measurement of in- 
terstellar dust occurs during the cruise phase as shown in flgure |^. During the 
short time of the encounter with the comet cometary dust particles are collected 
which are lifted from the comet's nucleus and accelerated into the coma by the 
drag of the outflowing gas. The amount of material in a comet's coma depends 
on the activity of its nucleus, which, in turn, depends on the distance from 
the Sun and on the properties of the nucleus' surface. To determine P/Wild 
2's activity we use ground based observations. For the properties of the dust 
phase, especially the mass distribution of dust particles, contained in a comet's 
coma, we rely on data collected with the VEGA and Giotto spacecraft at comet 
P/Halley during its visit of the inner Solar System in 1986. Since P/Halley and 
P/Wild 2 differ in activity and size, and since the VEGA and Giotto measure- 
ments had a different geometry than the measurements of Stardust, we have 
to construct a model of the coma and use the measured activity, brightness 
and mass distribution as an input to this model. By using the model to re- 
produce the spacecraft measurements of the dust fluence at P/Halley, we can 
validate the model and the procedure of taking the input activity and brightness 
from ground based observations. The secondary goal of the Stardust mission is 
the collection and in-situ measurement of contemporary interstellar dust parti- 
cles which enter the Solar System from the local interstellar cloud (LIC). The 
same instruments are used for the interstellar dust measurements as for the 
cometary part of the mission. During deflnite parts of Stardust's cruise, the 
dust collector will expose its backside to the interstellar dust stream and on 
other occasions CIDA will be pointed into the upstream-direction (see figure 
The interstellar dust stream was discovered by Griin et al. (1993) |Q in the 
data collected by the dust-detector on-board the Ulysses spacecraft after fly-by 
of Jupiter. It was clearly identified and distinguished from interplanetary dust 
by its opposite impact direction, its impact velocity in exceess of the escape 
velocity at Jupiter distance, and its constant impact rate at high heliocentric 
latitudes. The Ulysses measurements have been confirmed by impacts from 
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Figure 1: The trajectory of Stardust (solid line) in the ecliptic plane (the 
verenal equinox direction is to the right). The orbits of P/Wild 2 and Earth 
are indicated as dashed lines. Mission events are marked with the time given in 
days after launch in brackets. The flow direction of interstellar particles from 
the Lie as derived from the Ulysses/ Galileo measurements is shown by the 
arrow in the lower left corner. Thick lines indicate the part of the orbit used to 
collect interstellar dust with the Aerogel collector. 

the retrograde direction measured by the identical dust-detector on-board the 
Galileo spacecraft |^ . The total flux of interstellar particles was determined to 
be 1.5 • 10--* m-2 The mass distributions of the Ulysses and Galileo 

measurements have been found to be nearly identical and range from 10""'^^ g 
to 10~* g with particles of masses of 10~^^ g being most abundant Q. From 
the impact direction it was derived that interstellar dust particles come from 
the same direction Q as interstellar neutral Helium which was measured by the 
Ulysses/GAS experiment [Q. The upstream-direction of the interstellar He- 
lium was determined to be Aoo,gas = 254.7° ±1.3° (hehocentric longitude) and 
/3oo,gas = 4.6° ±0.7 (heliocentric latitude) |Q. Furthermore the relative velocity 
of the gas was given by Woo.gas = 25.4 ± 0.5 km s~^. These parameters relate the 
interstellar material measured in the Solar System with the LIC which has been 
identified in the Doppler shift measured in the UV by HST-GHRS The 
Doppler-shift indicates a relative motion between the Sun and the LIC with a 
velocity of 25.7 ± 0.5 km s^^. 

In this paper we use the expressions "flux" for the number of impacting 
particles of a given mass per unit area and time, "fluence" for the flux integrated 
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Figure 2: Equilibrium surface temperature versus the angle from the sub- 
solar point as given by the solution of equation (|^). 

over time, "total flux" for the flux of particles of all masses, and "total fluence" 
for the total flux integrated over time. 

2 Coma Model 

In the following we describe the coma model we use to determine the cometary 
dust flux on the Stardust spacecraft. Because there is no standard model of 
the gas and dust environment of a comet nucleus, yet, we describe the assump- 
tions and calculations we use in detail. So the reader can see how different 
assumptions or values of parameters change the results. 

Since the dust is dragged into the coma by outflowing gas, we first have 
to calculate the gas density and velocity inside the coma. We assume that 
the outflow is symmetric about the Sun-nucleus axis. Therefore, all spatial 
distributions depend only on the angle i!)s between the considered location and 
the axis of symmetry and the radial distance r. We furthermore assume that 
the dust particles do not affect the gas flow, i.e. we use a test-particle approach. 
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Table 1: Parameters for computation of the surface temperature. 

Solar constant /sun = 1370 W m^'^ 

nucleus surface albedo^ Ab ~ 6% 

Emissivity of comet surface e = 1 

Boltzmann constant k — 1.38 • 10~^^ J K 

Stefan Boltzmann constant a — 5.67 • 10~^ W K^'' m^'^ 

Advogadro Number Na — 6 • lO^^moP^ 

Latent heat of water ice L = 36 kJ mol~^ 

Vapor pressure parameter A = 3.56 • 10^^ N m~^ 

Vapor pressure parameter B — 6141 K 

^ bolometric Bond albedo 



2.1 Gas Phase 

To calculate the state of the gas around the nucleus, we need to know the activ- 
ity distribution on the nucleus' surface, which is governed by the temperature 
distribution. 

The temperature at a location on the day side of the nucleus' surface is 
calculated such that there is an equilibrium between solar illumination, water 
sublimation, and blackbody reradiation. Therefore, the surface temperature 
Ts{i}s) is the solution of the following equation (the meaning and values of the 
parameters in the next two equations are given in table hi) . 



i^'AB){^^^y -cosies) 



where ZhkI^s) is the Hertz-Knudsen sublimation rate given by: 



(1) 



ZMTs) = "^""'^ (2) 

We use the approximation p(rs) ~ Ae^^^^'^'^^ for the vapor pressure |lj]. On the 
night side of the nucleus the temperature is constantly set to 150 K. This value 
is also used at the day side near the terminator if the solution of equation 
is smaller than 150 K. Figure |^ shows the equilibrium temperature distribution 
determined from equation (0). 

Now, we obtain the activity Qhk of the comet by inserting the temperature 
distribution into equation and integrate over 'd^- 

Qhk = 2TrrlJ Z^Ki'd s)d cos {'ds)- (3) 

In our model we use the outgasing activity of the nucleus given by Farnham and 
Schleicher (1997) On 1997 March 5 they report an OH-activity Qoh = 5.9- 
10^^ s~^ and a value of the albedo-filling factor-radius product of Afp = 4.27 m. 
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Figure 3: Solid lines represent lines of constant gas density (arbitrary units) 
in the region up to 25 comet radii (= 50 km) around the nucleus. The Sim 
direction is upward and the short radial lines indicate the gas velocity (the scale 
on the top can be used to determine the modulus of the velocity). The gas 
density above the comet night side (bottom) is much lower than the density 
above the day side. However, the gas velocity above the day and night side is 
comparable. Apart from a region above the terminator, the gas flow is almost 
radial close to the nucleus. 



At this time the comet was at a heliocentric distance of 1.7 AU and was observed 
at a phase angle of 27°. Although the observation was taken on the in-bound 
part of the comet's orbit, and the Stardust spacecraft will approach the comet 
on the out-bound part, we assume that the observed activity is comparable to 
the comet's activity during the encounter. Therefore, we run our coma model 
for a comet at 1.7 AU with a water activity equal to the observed OH- activity 
QH2O = 5.9 • 10-^^ s~^. To match the observed activity, the water sublimation 
rate is calculated by scaling the Hertz-Knudsen sublimation rate: 

Zn.oi^s) = ^-ZukWs)) (4) 

In addition to the water activity, we assume a CO-activity of 10% of the wa- 
ter activity, i.e. Qco = 5.9 • 10^^ s~^, and a constant distribution above 
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the nucleus' surface Zcoi^s) — Qco/(47rr^). For calculating the accelera- 
tion of the dust particles by the cometary gas, we have to consider the most 
abundant gas molecules only, therefore other constituents than H2O and CO 
are not considered in the model and the total gas activity is set to Qgas = 
QH2O + Qco = 6.5 • 10^'' s~^. Accordingly, we set the total gas activity dis- 
tribution to Zs{ds) = Zii.^o{''^s) + ^co(i^s)- Using Zsi'&s), the thermodynamic 
state of the gas can be determined by a numerical procedure described by Crifo 
et al. (1995) § and KnoUenberg (1993) ||. 

As an input to the calculation of the gas state we use a mean mass of a gas 
molecule of mgas = ('TiHaoQHaO + wcoQco) /(QhsO + Qco) and a mean degree 
of freedom per molecule of /gas = (/H20QH2O + /coQco)/(Qh20 + Qco)- This 
simplification introduces some physical inconsistency into the model, but the 
results do not depend significantly on these choices. 

Another important parameter for modeling the gas outflow is the initial 
gas temperature above the nucleus' surface, because it determines the velocity 
reached by the gas, and therefore also has an impact on the dust acceleration. 
The initial gas temperature above the nucleus' surface does not necessarily co- 
incide with the surface temperature, because the gas is not at rest with respect 
to the surface. This effect was taken into account as it is described in Crifo and 
Rodionov (1997) || (Appendix B). 

As shown in figure ^ the gas density determined by our thermodynamic 
model is distributed highly asymmetric around the nucleus. On the day side 
much more gas is available to accelerate dust from the surface. The flow direc- 
tion is nearly radial. In the outer region of the coma, the gas is not dense enough 
to exert a significant drag force on the dust particles. We assume that the gas 
drag is negligible outside a maximum distance rmax of 25 times the radius of the 
nucleus. After leaving this region, the dust particles move independently from 
the gas. 

2.2 Dust Phase 

The dust phase of a comet's coma consists of all dust particles lifted and acceler- 
ated from the nucleus by the gas drag. The mass range of dust particles emitted 
by a comet covers many orders of magnitude and the dynamics of these particles 
depend on their mass. Therefore we divide the mass range [10~^° kg, 10° kg] 
into 20 logarithmic mass decades [rrii-i, rrii], i = 1, . . . , 20 which we call "dust 
classes" in the following. All dust particles which are released with a mass 
inside one of the intervals i are modeled by a dust particle of a representa- 
tive mass md,i. We assume that no processes which change the properties of 
a dust particle, e.g. fractionation due to mutual collisions, take place after a 
particle has left the nucleus' surface. Furthermore, we assume the dust parti- 
cles to be spherical with a constant bulk density = 1 g cm^'' and a radius 
Si — (3r7id^i/(47rpc;))'^^/'^^ For the determination of the amount of dust in the 
coma, we derive a dust-to-gas mass ratio Xi for each dust class. Using the dust 
to gas ratios Xi i we calculate the number of dust particles released per unit area 
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and time. 



ZdA^s) = X.^^gas(l?.) (5) 

md,i 

We take into account that particles cannot be lifted from the comet surface 
above a given mass by setting the dust activity to zero at locations on the 
comet surface where the gravitational attraction of the comet nucleus exceeds 
the drag force exerted by the gas. 

To calculate the gravitational force of the nucleus, we determined the nu- 
cleus' size by using ground-based observations of P/Wild 2 at large heliocentric 
distances. Fitzsimmons and Cartright (1995) (l^ report a red magnitude of 
i? = 22 at the heliocentric distance of r/j = 4.4 AU, a geocentric distance of 
A = 4.0 AU, and a phase angle of a = 12°. We assume that P/Wild 2's nu- 
cleus has a spherical shape. The radius of the nucleus r„ can be calculated by 
(compare Jewitt (1991) pef) 



A / l0-o-4(«-flsu„) 

1 AU y Pnuc • Jimc(a) 

where i?sun = —27.22 is the red magnitude of the Sun M, p^uc = 4% is used for 
the geometric albedo of the nucleus and jnuc(Q;) — 10"''^" is the phase function 
of the nucleus with S = 0.035 deg~^ |Q. We find r„ — 2.3 km. Because 
Fitzsimmons and Cartright (1995) [|6| reported a "near-stellar appearance", 
there was already some contribution of the dust coma to the observed intensity. 
Therefore, r„ = 2.3 km is an overestimation of the radius of the nucleus and we 
use the value r„ = 2 km in our model. 

Assuming a bulk density of /?„ = 0.5 g cm^"^, the gravitational acceleration 
is given by: 



Mn— T with /i„ = 1118m^s ^. (7) 
\r\ 



The acceleration of a dust particle due to gas drag can be expressed using the 
drag coefficient Cd- 

_ 1 ^lf]_r I- I 

^diag — (-j^gas^^gas ' l^gas ^dust| 

■ (Wgas - t/dust) (8) 

The drag coefficient in the free molecular approximation is given by Probstein 
(1968) M 



Cd - ^J^ + ^e(-^^) 



3s y Tgas ' S^y/lT 

45"* + 4s2 - 1 

- \^ erf(.), (9) 
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Figure 4: Trajectories of dust particles with radii of 8.8 /im, which start at 
equidistant positions on the surface. The scales in the image are nucleus radii. 
The velocity of the particles is dominated by the radial component, far from 
the nucleus the tangential velocity component is negligible. The distance of the 
trajectories above the terminator is increased due to the tangential gas flow in 
this region. 



where s = \vgas — Vdust\/ \/2fcTgas/?7T-gas- The dust temperature is set to Td = 
228 K (compare Divine (1981) 0: Ta = 310 K • (r^/l AU)-°-58 and rn = 
1.7 AU). At every point on the nucleus' surface for which the gas drag ac- 
celeration exceeds the gravitational acceleration of the nucleus, the equation of 
motion of a dust particle can be solved numerically inside the gas flow (see figure 
^ and For large particles which can be lifted only from the comet's day side, 
there are trajectories for which the radial velocity is not directed away from 
the nucleus but fall back into the nucleus' direction (see figure H). Since this 
behavior is only found in a narrow region, we do not consider these particles in 
the following. 

Apart from a narrow region near back- falling particles, the dust velocity 
is dominated by the radial velocity component. Since the contribution of the 
tangential velocity to the total dust velocity decreases with increasing distance 
from the nucleus' surface, the dust trajectories can be approximated by straight 
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Figure 5: Trajectories of dust particles with a radius of 8.6 mm. These big 
dust particles can only be lifted in a region around the; sub-solar point. Near the 
boundary of this region there arc trajectories which start to fall back in nucleus 
direction. In our model we neglect the contribution of these dust particles to 
the dust number density of the coma. 



radial lines with constant radial velocity at large distances. Outside the maxi- 
mum distance rmax we approximate the dust trajectories by straight, radial lines 
and constant velocities Ve^ii'ff) which is equal to the radial velocity computed 
for a dust particle at the distance Tmax- The angle i? denotes the angle of the 
position of the dust particle with respect to the Sun direction at rmax, and is 
also the angle of the extrapolated trajectory with the Sun direction. 

The dust number density outside the radius rmax can then be computed using 
the position of the trajectory at the nucleus' surface as a function of the angle 
■!? with respect to the Sun direction. To calculate the number density nd,i at 
positions with r > rmax and angle •&, we introduce the activity distribution of 
the comet outside the acceleration zone: 



^dcos{'&s) 
dcos{'d) 



Z,, = rlZM^))^^^^^ (10) 
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For the number density n^.i we get: 



The number density inside the acceleration cone can be calculated in an analo- 
gous way. 

As a result, we show the contours of constant dust particle density in figure 
|. The region of high dust density occupies a larger volume on the day side 
than on the night side. Due to the drop in activity at the terminator, the dust 
particle density right above the terminator is lower than on the night side. Dust 
particles from the day side can propagate freely to the night side which extends 
the region of high dust particle density just behind the terminator in a narrow 
radial feature. If this is a real feature at a comet and if it can be observed in 
reality remains an open question. The dust particle density opposite to the Sun 
direction is the result of the night side activity. 

We want to use the model to determine the dust flux on the Stardust instru- 
ments during the fly-by of the nucleus. We can neglect the influence of radiation 
pressure on the particle's trajectories, because Stardust approaches the nucleus 
sufficiently close. Because the velocity of Stardust relative to nucleus is much 
larger than the velocity of the dust particles, the dust impact velocity at the 
spacecraft is dominated by the spacecraft's motion. We therefore calculate the 
dust flux on the spacecraft by the product of the dust number density and the 
spacecraft velocity wsc relative to the nucleus. 

/sc = na^, ■ vsc (12) 



From the result in equation ( |12D we can derive the dust fluence on the instru- 
ments, which is given by the time integral over the flux and represents the 
column density of dust particles along the spacecraft's trajectory. We need the 
column density of dust particles also to calculate the brightness measured by 
an ground-based observer for a given line of sight. 

As we show in Appendix A, the column density (in number of dust particles 
per m^) for an image plane with phase angle a can be calculated by 

— ' dcos(v) 



d{x',y') pj_^ VeA"^) 

•^=^===A^(a,r;,^), (13) 

y'U'' — (COS(W)^) 



where u = u{a, rj) = -^/l — (sin(r7) • sin(a;))^. The variables p and ry represent 
the position of the line of sight in the projected plane: p is the distance from 
the (projected) nucleus position and ry denotes the direction of the hue of sight 
position with anti-solar direction projected in the image plane (compare figure 
pO|) . N(a, rj, -d) is equal to the number of intersections of a line of sight with the 
cone of constant i!) and is equal to 0, 1, or 2 (see table ||). 
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Figure 6: The lines of constant dust number density of particles with radii of 
&.8iim in the region within 25 radii of the nucleus. The radial lines indicate the 
velocity of the dust particles. Because of the low gas density, the dust particle 
velocity is smaller above the night side. Due to the increased tangential distance 
of the dust trajectories, the dust number density is small above the terminator 
(see figure 0). 



For comparison with our axially symmetric model we provide the column 
density in the radially symmetric case: 



dUrl 



d{x',y') 



PVe,i 



(14) 



2.3 Adjustment of the Dust Activity to the Observed Afp- 
value 

We adjust the dust activity of our model so that the observed Afpohs = 4.27 m 
is matched. The ^/p-parameter is related to the intensity /dust received from 
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Table 2: Number of intersections of a line of sight of angle rj with a cone 
of opening angle ■& for different cases as explained in the text. Note that 
(A^(a,r?,i?) + Af(a,7r- r/,i9))/2 1 for |cos(i9)| < u. 



Condition 1 


Condition 2 




u < cos('i9) 


none 





cos(a) < cos(i9) < u 


cos(?7) < 
cos(?7) > 


2 



cos(a) > cos{-&) > 
— cos(a) 


none 


1 


— |cos(a)| > cos(i9) > —u 


cos{rj) < 
cos(?7) > 



2 


—u > cos{-d) 


none 






dust particles inside a field of view with radius p by A'Hearn et al. (1984) 



Afp . 4^!%^.^, (15) 



Sun 



where /gu„ is the radiation flux density of the Sun at a heliocentric distance 
of 1 AU observed in the same band of the spectrum as the comet, A is the 
geocentric distance and r/j the heliocentric distance. 

The intensity received due to a single dust particle is given by 



where p is the dust geometric albedo and j{a) is the dust phase function of the 
phase angle a. Using the number of dust particles inside the aperture Nd.i (see 
equation (|3^)), we can determine the ^/p-parameter of our model by 



/'27r /'U 

^/Pmod = 4pj(a)Vs? / / 

Jo J ~u 



yjv? — C0s2(l?) 



dcos{^)dr]. (17) 



Here we assume that all dust particles have a geometric albedo and a phase 
function independent of size. For the geometric albedo we assume a value of 4% 
and take the phase function from Divine (1981) pO[ . 

To match the observed Afpohs with the coma model, we have to scale the 
dust activity of the comet, i.e. the dust to gas ratio of all dust classes have to 
be scaled by the factor A/pobs/^/Pmod- 

2.4 Dust Mass Distribution 

The mass distribution of the dust particles is the most crucial parameter of 
the model, because a change in the mass distribution can change the total dust 
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fluence by orders of magnitude. Therefore we take special care to choose a 
reasonable mass distribution. In this section we explain our choice and discuss 
the resulting fluence in comparison with the measured fluence in the coma of 
P /Halley. The reader might want to check what the effects of a different choice 
of the mass distribution are. 

The VEGA spacecraft measured the number of particles per unit area as 
a function of mass. We use this measurement as the mass distribution of our 
model. As an analytical representation of the mass distribution, we use a fit by 
Divine and Newburn (1987) 



The parameters which fit the cumulative mass distribution of the VEGA 2 
fluence best are a = 0.9, b = 0.29, c = 2.16, and to* = 1.6 • lO^^^ kg 0. Since 
the mass distribution is only used to determine the relative but not the total 
abundance of particles in different dust classes. Ft can be set to an arbitrary 
value. Note that —a — —0.9 is the exponent of the cumulative mass distribution 
for big dust particles (m ^ mt,x ^ 1), and — ac — —1.9 is the exponent for 
small particles (m <C mt, x ^ 1). 

For each dust class [mi_i,mi], i — 1, ... ,20 considered in the model, the 
fluence on the VEGA 2 spacecraft inside one of the mass bins is proportional to 



The dust fluence on the VEGA 2 spacecraft can be modeled analogous to 
the fluence on Stardust. At closest approach, VEGA 2 was also far inside 
the envelope of the measured dust particles, therefore, the total fluence was 
dominated by direct particles. 

Because of the dependency of the integrand in equation (|3|) on the dust 
escape velocity Ug^i, which, in turn, depends on the mass of the dust particle, 
the mass distribution of particles which hit the VEGA 2 instruments does not 
coincide with the mass distribution of particles at the nucleus' surface, i.e. the 
distribution of Zd^i- Therefore a correction has to be applied to the VEGA 2 
fluence before it can be used as mass distribution at the nucleus' surface. For 
this purpose we compute the dust escape velocities We,i,HaiiGy for an analogous 
radially symmetric P/Hallcy model (The model parameters are QH2O — 6.8 • 
^q29 g-i ^ggg Krankowsky (1986) @, Qco — s~^, r„ = 5 km; all other model 
parameters are not changed). The fluence of dust class i is inversely proportional 
to the dust escape velocity Ve^i_HaUey for a radially symmetric coma model (see 
equation (|^). Therefore, the product Usr/.i = 7^^/G2,^ ■ Ve,i,HaUoy is proportional 
to the mass distribution at the comet's surface and the dust-to-gas mass ratios 
of the dust classes can be computed by 




(18) 



nvG2,i oc F{mi^i) - F{mi) 



(19) 



X2 — Xaux 



(20) 
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where Xaux is an auxiliary dust to gas ratio. Xaux would have the meaning of the 
ratio of dust mass to gas mass released per unit time from the comet's nucleus, 
if all dust classes could be lifted from the comet's surface. Since there are dust 
classes which can not be lifted from the comet's surface or can only be lifted 
inside some region around the sub-solar point, no physical interpretation can be 
linked to this number and is in fact only an auxiliary value needed during the 
computation. 

3 Model of Lie Dust Traversing the Solar Sys- 
tem 

The impact parameters of interstellar particles on the dust collector and the 
CIDA target are determined by the impact direction and velocity of the particles, 
and also by the articulation of the spacecraft during the cruise phase. Since 
both is non-trivial, we separate the effects by discussing two different types of 
particles: (1) the reference particles as defined by the Stardust mission plan |^ 
which move on straight lines through the Solar System in a direction derived 
from Ulysses and Galileo data, and (2) larger particles coming from the same 
direction, but being deflected due to solar gravity. For the reference particles 
the change in impact parameters is caused only by the spacecraft's trajectory 
and the articulation strategies defined in the mission plan for the cruise phase. 
Larger particles are mainly interesting for the dust-collector, because CIDA is 
unlikely to be hit by such a particle due to its small target area and the low 
abundance of these particles. Furthermore, larger particles can be more easily 
extracted from the aerogel, this extraction might be difficult for particles with 
sizes of about a few tenth of a micron. In the following we argue why the 
reference particles are a good approximation of the most abundant interstellar 
particles, which are to be measured with Stardust and why smaller interstellar 
particles are not expected to reach the spacecraft. 

In general, the dynamics of a charged, massive dust-particle in the inter- 
planetary environment can be described by 

f = -7(l-/3)^-h£(z/.wf xs). (21) 

This equation of motion takes into account gravity, radiation pressure, and 
the Lorentz force induced by the solar wind magnetic field. 7 is the coupling 
constant of gravity, Mq is the solar mass and the parameter f3 is defined as 
the ratio of the magnitudes of radiation pressure force and gravity. Q is the 
equilibrium charge to which a dust-particle will be charged in the interplanetary 
radiation and plasma environment, Wg^f is the velocity vector of the particle, 
measured in the frame of the radially expanding solar wind, and B represents the 
solar wind magnetic field. In equation (|l]) drag forces like Poynting-Robertson 
drag and solar wind drag are neglected, because these drag forces are small 
and only have a long-term effect on particles on bound orbits. The dynamical 
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parameters (5 and Q/m depend on the particle size and can be represented as 
functions of particle radius a, if we assume a spherical shape for the particles 



3.1 Reference Particles 

Physically, the reference particles are implemented by setting /? = 1 and Q jm = 
C kg^^. In this case the solution to equation ( [2l] ) are trajectories along straight 
lines, i.e. r — '?oo,dust = const.. Since the Ulysses measurements indicate that 
interstellar dust and gas are kinematically coupled, we set Woo.dust = "^cxj.gas using 
the parameters derived for the gas by Witte et al. (1996) Q Gustafson (1994) 
[ p3| calculated /3 as a function of particle radius assuming a spherical shape and 
a composition like the astronomical silicates defined by Draine and Lee (1984) 
[ p!2| (see figure]^). Using these calculations it was shown Q, that for 57% of 
all interstellar particles measured by Ulysses and Galileo |1 — /3| < 0.4. If the 
particles are measured at about 2 AU, this deviation from the reference behavior 
would change the impact angle by less than 15°, which is the uncertainty in the 
determination of the interstellar dust upstream direction from the Ulysses and 
Galileo measurements. The impact velocity will deviate by less than 30% for 
these particles. Furthermore, the variation of impact direction and velocity 
due to electro-magnetic interaction of interstellar particles with the solar wind 
magnetic field [Q/ra 7^ C kg""'^) can be neglected if they have radii larger than 
0.2 /im. 



1.000 



0.100- 



0.010 r 



0.001 
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Figure 7: The parameter /3 as a function of particle radius a. Homogeneous, 
spherical particles with the optical constants of astronomical silicates have been 



assumed (see [23 ) 
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3.2 Gravitational Particles 



As shown in figure « 1 is not valid for large particles. Large interstellar 
particles (m > 10~^^ g) are less abundant than smaller (m < 10~^^ g) ones 
[H, so they are more important for the dust collector than for CIDA due to 
its larger surface area. If we neglect electro-magnetic effects, the trajectories 
are hyperbolae and the spatial density distribution of these particles can be 
calculated analytically (see equation ( p] ) in Appendix A). Figure ^ shows this 
distribution as a grey-scale in the ecliptic plane. In this simple model the density 
enhancement of large particles due to gravitational focusing during the collection 
phase is less than a factor of 2. The enhancement in spatial density downstream 



4 




AU 

Figure 8: The normalized density distribution of gravitationally deflected 
(/? — 0.1) interstellar particles in the ecliptic plane represented by a grey-scale. 
The normalization is defined to set the spatial density at infinity equal to 1. 
Stardust's trajectory is indicated as the white solid line (compare to figure 1). 

the Sun is caused by gravitational focusing. Since the downstream vector of the 
interstellar dust flux has an angle of « 5° with respect to the ecliptic plane, 
the actual focus does not lie in the ecliptic, but a peripheral spatial density 
enhancement is present in the ecliptic plane. The main effect of the dynamics of 
large particles will be a higher impact velocity and a deviation in impact angle 
with respect to the reference particles. 

Particles with large /3- values (/3 > 1.4) have diameters of w 0.45 /^m which 
correspond to the maximum wavelength of the solar spectrum. As we will argue 
in the next section, we do not expect these particles to be abundant in the Solar 
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System during the Stardust mission due to their interaction with the solar wind 
magnetic field. 



3.3 Electromagnetic Effects 



In section B.l we have argued that electromagnetic effects can be neglected for 
the determination of the impact direction and velocity of interstellar particles. 
But we have to take into account these effects when considering the abundance 
of interstellar particles in the Solar System since long-term effects might reduce 
or enhance their spatial density. It was argued by Levy and Jokipii (1976) 
[ p6| that classical (very small) interstellar particles are removed from the Solar 
System by Lorentz force. The Ulysses and Galileo measurements show |^ a de- 
pletion of small particles (but still one order of magnitude above the detection 
threshold). Models of the electromagnetic interaction of the charged interstel- 
lar particles with the solar wind magnetic field ||2l|, ^ |l9| predict a periodic 
focusing and defocusing of the particles to the solar equator plane with the 22- 
years solar cycle. A new magnetic cycle starts at solar maximum and the mean 
deflection effect is strongest when the magnetic field is in the most ordered con- 
figuration during the solar minimum. The last solar maximum in 1991 started 
a defocusing cycle. We have calculated the total flux of interstellar particles on 
the Ulysses detector as a function of time to check if the defocusing causes the 
total flux to drop. Figure ^ shows that the total flux (in the heliocentric iner- 
tial frame) of interstellar particles was constant after Ulysses' fly-by of Jupiter 
in February 1992 (when Ulysses left the ecliptic plane). The gap in the data 
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Figure 9: Total flux / of interstellar particles measured by Ulysses as a function 
of time. The error-bars are due to the fraction of particles which could not be 
clearly identified as interstellar dust. 



between beginning and mid of 1995 is caused by the ecliptic crossing in March 
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1995 when distinction between interstellar and interplanetary particles was dif- 
ficult. After mid of 1996 the measured total flux of interstellar particles drops 
by a factor of 3. We interpret this phenomenon as the effect of electromagnetic 
defocusing. Unfortunately this can not be confirmed by Galileo data, because 
Galileo is inside the Jovian system since the end of 1995, therefore, Jupiter dust 
dominates the data and interstellar dust impacts can not be identified |Q . We 
assume that the period of reduced total flux will last for half a solar cycle, which 
is 11 years, i.e. from mid of 1996 to mid of 2007. This period of time covers 
the duration of the Stardust mission. Following this estimate, we predict from 
the Ulysses data that the total dust flux on the Stardust experiments will be 
reduced by a factor of 3 compared to the Ulysses and Galileo measurements. 

For Ulysses and Gahleo a total flux of /u/g — 1-5 • 10^"* m^^ s^^ has been 
determined so we predict /sd = 5 • 10~^ m^^ s^^. 

4 Results 

4.1 Cometary Dust Measurements 

In this section we apply the coma model developed in section ^ to quantitatively 
predict the flux and the fluence of cometary dust particles on the Stardust 
spacecraft during the P/Wild 2 encounter. The results we obtain for an axially 
symmetric model are compared with the results we get if we use a radially 
symmetric model. We can compare our results for Stardust at P/Wild 2 with 
the fluence measured by the VEGA spacecraft inside the coma of P/Hallcy 
by scaling the predicted fluence with the difference in brightness, encounter 
distance, and phase angle. 

Fluence as a function of mass We approximate the trajectory of Stardust 
as a straight line passing the coma with a phase angle of 70° with respect to the 
Sun direction. Stardust's closest approach to the nucleus is above the comet's 
day side. The fluence of dust particles of a given dust class i during the whole 
trajectory is equal to the column density along the spacecraft's trajectory which 
is given by equation ( p^ ) . As the fluence is inversely proportional to the distance 
of closest approach of the spacecraft to the nucleus p, our results can easily be 
scaled to other closest approach distances. The numerical values of the fluence 
for a trajectory above the sub-solar point is given in table |. 

Flux and Fluence as a lYinction of Spacecraft Position The fluence as 
a function of time from closest approach is calculated by the time integral, or, 
equivalently, by the integral over the path s along the spacecraft's trajectory, of 
the dust flux on the spacecraft as given in equation ([l^). 




(22) 
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Tabic 3: Flucncc of dust particles on the spacecraft for the radial and axially 
symmetric model for a closest approach distance of p = 100 km. For different 
distances p' of closest approach the fluence can be scaled hy p/ p' . In the axially 
symmetric model we assume that the spacecraft's trajectory crosses the sub- 
solar point. 



Dust 


particle 


particle 












class 


mass 


radius 






fluence 


[m-2] 




i 


md,i [kg] 


Si [m 






radial-symm. 


axial- 


symm. 


1 


10-20 


1.34- 


10" 


-8 


5.86- 


10+uv 


1.11- 


10+08 


2 


10-19 


2.88- 


10- 


-8 


3.86- 


10+07 


7.05- 


10+07 


3 


10-18 


6.20- 


10- 


-8 


2.52- 


10+07 


4.46- 


10+07 


4 


10-17 


1.34- 


10- 


-7 


1.60- 


10+07 


2.76- 


10+07 


5 


10-16 


2.88- 


10- 


-7 


9.67- 


10+06 


1.64- 


10+07 


6 


10-15 


6.20- 


10- 


-7 


5.39- 


10+06 


8.98- 


10+06 


7 


10-14 


1.34- 


10- 


-6 


2.58- 


10+06 


4.26- 


10+06 


8 


10-13 


2.88- 


10- 


-6 


9.37- 


10+05 


1.53- 


10+06 


9 


10-12 


6.20- 


10- 


-6 


2.32- 


10+05 


3.75- 


10+05 


10 


10-11 


1.34- 


10- 


-5 


4.11 • 


10+04 


6.62- 


10+04 


11 


10-10 


2.88- 


10- 


'5 


6.02- 


10+03 


9.64- 


10+03 


12 


10-09 


6.20- 


10- 


-5 


8.11 • 


10+02 


1.29- 


10+03 


13 


10-08 


1.34- 


10- 


-4 


1.06- 


10+02 


1.67- 


10+02 


14 


10-07 


2.88- 


10- 


'4 


1.37- 


10+01 


2.14- 


10+01 


15 


10-06 


6.20- 


10- 


-4 


1.78- 


10+00 


2.72- 


10+00 


16 


10-05 


1.34- 


10- 


-3 


2.40- 


10-01 


3.44- 


10-01 


17 


10-04 


2.88- 


10- 


-3 


3.55- 


10-02 


4.34- 


10-02 


18 


10-03 


6.20- 


10- 


-3 


6.31- 


10-03 


5.37- 


10-03 


19 


10-02 


1.34- 


10- 


-2 


0.00- 


10+00 


5.71- 


10-04 


20 


10-01 


2.88- 


10- 


-2 


0.00- 


10+00 


1.15- 


10-05 
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For a radially symmetric model an analytical expression is given by: 



F. 



SC.i 



Vsct 



ds 



(23) 



where ^ = arctan(VscVp) is the angle of the spacecraft position to the point of 
closest approach. From equation (|2^) we see that Fsc,i is a linear function of ^ 
in the radially symmetric model. 

We show the calculated total dust fluence on Stardust in figure nu. In the 



time from closest approach 

-Imin -Ijs -^5 Q 15s ln;iLn 
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Figure 10: Total fluence on Stardust versus the angle ^ from the location of 
closest approach for the radially symmetric (dashed line) and axially symmetric 
model (solid line) . The total fluence calculated by the axially symmetric model 
deviates significantly from the result of the radially symmetric calculation. 

axially symmetric model the total fluence is approximately a linear function of 
^ for —40° < £. < 40° , like it would be for the radially symmetric model. The 
radially symmetric model predicts a too low dust flux during closest approach, 
because it does not take into account the enhanced activity on the day side. 
Therefore, the slope of the total fluence is flatter between ^ — —60° and ^ = 60° 
in the radially symmetric case. For ^ > 60° Stardust enters the region above the 
terminator and consequently the total flux (and therefore the slope of the total 
fluence) drops below the value predicted by the radially symmetric model. We 
conclude that, for predicting the total dust fluence on Stardust during closest 
approach, one must take into account that Stardust approaches the comet from 
the day side and passes the terminator shortly after closest approach. The 
different results from the radially and axially symmmetric calculations show 
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Figure 11: Total dust flux versus the angle ^ between Stardust's location on 
the trajectory and the location of closest approach P/Wild 2. The total flux of 
the axis and radially symmetric model are represented by the solid and dashed 
lines, respectively. It can be seen that the total flux predicted by the radially 
symmetric model is symmetric with respect to closest approach and reaches its 
maximum at that point. 



that the total fluence is strongly affected by the activity distribution on the 
surface of the nucleus. Unfortunately, we do not have any information about the 
activity distribution on P/Wild 2, so we have to rely on the axially symmetric 
coma model as a simple yet complete physical model of the inner coma for which 
we can constrain the parameters by measurements. 

In figure |ll| we show the total particle flux versus f . In the radially symmetric 
model the closest approach coincides with the point of highest total dust flux. 
Using the axially symmetric model, we find that the maximum total flux is 
reached 3 s before closest approach and that the total flux stays nearly constant 
until after closest approach. This is because the spacecraft enters into a region 
of lower activity after it has passed the sub-solar point but the decrease of 
cometary activity is apparently compensated by the still decreasing distance 
from the nucleus. 

Reproduction of the VEGA and Giotto measurements at P/Halley 

To validate our prediction of the fluence on the Stardust instruments, we scale 
our model to the conditions at P/Halley during the fly by of the Giotto and 
VEGA spacecraft. We stress that, although the VEGA 2 data were used to 
derive a mass distribution of our coma model, the comparison of the fluence 
is in fact a check of the model. This is because the VEGA 2 data were only 
used to derive the relative abundance of particles in different dust classes, but 
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we derive the total amount of the dust activity from the value of Afp that has 
been determined from ground based observations. 

The fluence on the Stardust and VEGA spacecraft is proportional to the 
value of Afp of the coma, inversely proportional to the distance p of closest 
approach, and the dust phase function during the observation. Thus, the fluence 
which was computed for the Stardust spacecraft can be scaled to the fluence 
measured on board the VEGA spacecraft by multiplying the fluence with the 
following three ratios: 

closest approach distance ratio: 

Pstarduat ^ 100 kill _ J_ 

PvBGA 8000 km 80 

A/p-ratio: 

Afpn^iioy ^ 219 m ^ gj^ o (24) 

AfPwild 1 ^ 4.27 m ~ 

dust phase function ratio: 
j(Q^wiid 2 ~ 30°) _ 0.04 



J (anaiicy ~ 60°) 032 



= 1.25 



The Aj p value of P/Halley during the VEGA encounter was taken from Schle- 
icher et al. (1998) Thisleadstoatotalscalingfactorof 51.3- 1.25/80 « 0.8. 
Hence the fluence which we predict for Stardust at P /Wild 2 is compatible with 
the fluence measured on the VEGA spacecraft at P /Halley. In figure ^ we show 
the data of the fluence measured by VEGA together with the our prediction of 
the fluence on Stardust scaled by 0.8. The good agreement of the number of 
dust particles collected on the VEGA spacecraft with the scaled Stardust flu- 
ence validates our procedure to derive the total dust activity of the comet from 
the observations and it shows that the parameter set we use in our model is 
consistent. 

Unfortunately, because of the loss of the Giotto data near closest approach 
due to large particle impacts, we can not perform any quantitative comparison 
of the Giotto fluence to our model. However, we can estimate the total mass 
fluence on Giotto on the basis of the well known deceleration of Ay — 23 cm s~^ 
of Giotto during the encounter. As for the considered mass distribution, the 
total fluence is dominated by large particles and the momentum enhancement 
factor is believed to be rather low for large particles Q , the total mass fluence 
is estimated by assuming inelastic dust particle impacts on Giotto: With the 
Giotto mass of Mciotto = 573.7 kg and the velocity of Vciotto = 68.37 km s~-^ 
during encounter, the total mass which hit Giotto can be estimated to be 
WGiotto = -^Giotto/ Vciotto ■ Ay « 2 g. Using the area of the Giotto dust shield 
oi A — 2.64: m^, we determine the total mass fluence on Giotto to be 0.76 g m~^. 
The total mass fluence on the Stardust at P/Wild 2 can also be scaled to the 
total mass fluence on Giotto at P /Halley by taking into account that the closest 
approach distance of Giotto was approximatively 600 km. We find the values 
0.038 g and 0.064 g for the radially symmetric and axially symmetric model, 
respectively. Thus, the extrapolation of the mass distribution which was fitted 
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Figure 12: Fluence of dust particles on the VEGA spacecraft. Filled and 
outlined triangles represent VEGA 1 and VEGA 2 data, respectively. The fit to 
the VEGA 2 data by Divine and Newburn (1987) ||ll[ is shown as the dashed 
line and the solid lines represent the scaled flucnces as predicted by the axially 
symmetric (upper solid curve) and radially symmetric model (lower solid curve) 

by Divine and Newburn (1987) |l^ to the fluence on the VEGA 2 spacecraft 
to larger masses leads to an underestimation of the total mass fluence by about 
one order of magnitude when compared to the total mass fluence derived from 
the final impact on Giotto. We discuss this discrepancy and its significance in 
section ||. 

4.2 Interstellar Dust Measurements 

We predict the measurements of interstellar dust during the interplanetary 
cruise phase by both main instruments on-board Stardust, CIDA and the aero- 
gel collector. Unlike in the case of the measurements at P/Wild 2, the geometry 
of the measurements of both instruments is different, because they are used at 
different parts of the interplanetary trajectory. Therefore, we discuss the pre- 
dictions for both instruments separately. 

CIDA To comprehend the predictions we make for CIDA one has to under- 
stand the three attitude strategies which have been defined by the mission plan 
[ p5| . Since CIDA has a fixed position on Stardust, its pointing is determined by 
the spacecraft attitude. The overall strategy should be to choose the attitude in 
a way that a maximum number of interstellar particles are detected assuming 
that they all behave like the reference particles. But the more one optimizes 
the attitude, the more complicated the spacecraft operations will become during 
the cruise phase. Thus, a trade-off has to be found between maximum number 
of detected particles and operation complexity. 
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We now explain the three different attitude strategies defined in the mission 
plan. Figure |l^ shows a simplified sketch of Stardust containing the high gain 
antenna defining the +z-direction , the solar panels, which lie in the x-y-plane 
and the CIDA field of view (FOV) in the x-z-plane. In this plane the angle 



Figure f3: Configuration of the CIDA field of view (FOV) with respect to 
Stardust. The +z-direction is defined by the pointing of the high gain antenna 
and the fixed solar panels, which should always be pointed roughly towards the 
Sun. 

of FOV is 52°. Strategy #1 fixes the +a;-axis to the upstream-direction of 
reference particles. This means that the +z axis, i.e. the solar panels and 
the high gain antenna, are not always pointing to the Sun. The mission plan 
restricts this off-pointing to 45° due to power concerns. If more than 45° off- 
pointing is required to keep the -l-x-axis pointing upstream, the spacecraft is to 
be turned to full Sun-pointing again, which means that the upstream-direction 
is turned out of the FOV. Strategy #2 is the most complex one. If pointing of 
the -|-a;-axis towards the upstream direction is possible without more than 45° 
off-pointing of the -|-z-axis with respect to the Sun, strategy #2 is identical to 
strategy #1. But if this configuration is not possible any more, strategy #2 
keeps the 45° off-pointing as long as the upstream-direction lies inside the FOV. 
If even this is impossible, the spacecraft should return to Sun-pointing again. 
The most simple strategy is strategy #3. It fixes Sun-pointing of the -f z-axis 
all the time, so CIDA will only collect interstellar particles when the upstream 
direction of the flux is occasionally in the FOV. Table ^ summarizes the three 
attitude strategies. 

Since the reference particles have a constant velocity of 26 km s~^, the 
impact velocity of the reference particles as a function of time is determined 
by the spacecraft velocity. The impact velocity, the impact energy, and the 
quantity mvf'^, which is proportional to the impact charge for impact plasma 
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Table 4: Attitude strategies for the Stardust mission during the cruise phase 



as given by the mission plan. 


strategy 


description 


# 1 


+a:-axis points upstream. Maximum 




Sun ofF-pointing is 45°. 


#2 


upstream-direction is kept in the FOV. 




Maximum Sun ofF-pointing is 45°. 


#3 


+z-axis always points towards the Sun. 



detectors p^ , where m is the mass and Vi is the impact velocity of the particle, 
are shown in figure 

We predict the impact rate v of reference particles on the CIDA target. We 
can estimate f by the given effective sensitive target area of ^cida — 80 cm^ • 



cos 40° = 60 cm^ the flux /sd = 5 • lO''^ m'^ s^ (see section in the 
inertial heliocentric frame, and the enhancement due to the upstream motion 
of Stardust. We get 

i^max - AciDA/sDf^^^^')«7.0.10-^S-\ (25) 

where Wrchmax = 60 km s^^ is the maximum relative velocity of Stardust to 
the stream of reference particles (see above) and Wdust,Gci = 26 km s^^ is the 
velocity of the reference particles in the inertial ecliptic frame. Figure |l^ shows 
the predicted impact rate on the CIDA target as a function of time for all three 
strategies explained above. 

The total number of reference particles that CIDA will detect is given by 
the accumulated impact rate over time. We show the prediction in figure |l6|. 

In summary, we predict that, assuming attitude strategy #2, CIDA will 
detect about 25 particles which behave very much like reference particles and hit 
the target with velocities between 30 km s^^ and 60 km s~^ and energies between 
18 eV and 4 eV N'^ (eV N'^ = electronvolt per nucleon) . For a summary 
of the prediction see table |5|, this table also contains the number of interplanetary 
dust particles as determined from the five-population model of interplanetary 
dust by Staubach et al. (1997) [E2| "contaminating" the measurement. 
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Table 5: Number of interstellar and interplanetary particles predicted to be 
measured by Stardust during the cruise. For the collector we give the numbers 
for the front (which is exposed to the comet) and the back side (which is exposed 
to the interstellar upstream direction) seperately. 
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^ Attitude strategy #2, see table 
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Figure 14: The values oi v, E — mv'^/2, and mv^'^ as a frmction of time, 
where m is the mass and v the velocity of the impacting reference particle. The 
shaded areas labeled "Arogel" indicate the phases of deployment of the Aerogel 
collector and "encounter Phases" labels the phase of preparations of the comet 
encounter. During these phases, CIDA is not taking data. 
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t [days after launch] 



Figure 15: Impact rate on the CIDA target as a function of time. Prediction 
in panel (a) assumes strategy #1, in panel (b) strategy #2, and in panel (c) 
strategy #3. Mission phases during which CIDA is not allowed to be turned on 
due to other activities are indicated by the shaded areas (compare figure |lj). 
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Figure 16: Total fluence on the CIDA target as a function of time. Prediction 
in panel (a) assumes strategy #1, in panel (b) strategy ^2, and in panel (c) 
strategy #3. Mission phases during which CIDA is not allowed to be turned on 
due to other activities are indicated by the shaded areas (compare figure 



29 



500 1000 1500 2000 2500 

t [days after launch] 



Figure 17: Total fluence on the aerogel detector of reference particles (a) and 
gravitationally deflected particles (b) as a function of time. Aerogel collection 
is not possible in the hatched periods. 

The Dust Collector For the dust collector larger particles are more in- 
teresting, because their trajectories can more easily be measured and they can 
be extracted from the aerogel. Since the collector is turned to maintain a point- 
ing to the upstream direction of the reference particles, we expect a spread in 
impact direction due to dynamical effects. For larger particles the dynamics 



should be dominated by gravity and we assume P — 0.1 (see section 3.2). The 



Mie-calculations |23| give a radius of 1 /im for these particles, which translates 
to a mass of 1 • 10^^^ g assuming a density of 2.5 g cm^"^. In the Ulysses and 
Galileo data particles heavier than 10"^"'^ g make up « 10% of the total number 



but they are not depleted due to electro-magnetic effects (see section 3.3) 



So we assume a flux of 1.5 • 10 ^ m ^ s ^ of these large particles. We show 
the total fluence of interstellar particles on the collector for both, reference and 
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500 1000 1500 2000 2500 

t [days after launch] 

Figure 18: Deviation of the impact angle of gravitationally deflected particles 
from the impact angle of reference particles, which should be close to 0° as a 
function of time. Aerogel collection is not possible in the hatched regions. 

gravitationally deflected particles, in figure |l^. 

Due to particle dynamics, tracks of larger particles should show a direction 
distribution. We show the deviation of impact angle as a function of time in 
figure 

For the collector it is furthermore interesting how fast gravitationally de- 
flected particles impact into the aerogel. Since the particles are accelerated 
towards the Sun we expect their impact velocities to be higher than the impact 
velocities of the reference particles. In figure |l^ we show the same plot as in 
figure ^4] for gravitationally deflected particles. 

To summarize our prediction for the aerogel collector, we give the total 
number particles collected during the 290 days of exposure time to be 120. 
80 of these particles are small, like reference particles. The remaining 40 are 
large particles, the impact tracks of which have angular deviations between 10° 
and 30°. For the large particles we expect impact velocities on the aerogel 
between 20 km s^^ and 40 km s~^. 
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Figure 19: Values v, E = and mv^'^ as a function of time, where m 

is the mass and v the velocity of the impacting gravitationally deflected particle 
(compare figure |lj) . 
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5 Discussion 



We think that the Stardust mission to comet P /Wild 2 will enhance our under- 
standing of the environment and the properties of comets by a large amount. 
According to our prediction, Stardust is able to fulfill its primary goal, the in 
situ measurement and collection of cometary dust particles. Concerning the 
secondary goal, the in situ detection and collection of interstellar dust particles 
from the local interstellar cloud, we find that the total number of particles we 
expect to be measured by Stardust is quite small considering the seven year 
duration of the mission. But having material from beyond the Solar System in 
a laboratory on the ground would be a big achievement. 

We have determined the fluence on the Stardust spacecraft using a coma 
model. In our model the dust activity was determined using observations of 
P/Wild 2 at 1.7 AU heliocentric distance. Our model can reproduce the mea- 
surements of the VEGA spacecraft at P /Hallcy when scaling the model to the 
geometry of this measurement. The total mass fluence on Giotto, deduced from 
the fact that the spacecraft was hit by a very large particle, is underestimated 
by our model by one order of magnitude. Of course, a one-particle hit has not 
a very high statistic significance and we can hypothesize that this one hit was 
a low probability event. The possibility that the mass distribution used in our 
model has a deficit of large dust particles, as suggested by the DIDSY data JTsf , 
exists, but simply adding large particles to the mass distribution would cither 
contradict ground based observations of the total brightness of P /Wild 2 during 
its in-bound part of its orbit at 1.7 AU, or the dust fluence measured by VEGA 
at P/Halley. The real mass distribution of dust particles at the nucleus is still 
unknown and the reader may choose a different mass distribution and put it 
into our calculation to calculate the expected fiux on Stardust. 

Stardust approaches the comet from the sunlit side and therefore more dust 
particle impacts are expected before than after closest approach. Thus, a ra- 
dially symmetric model is not sufficient to calculate the dust fluence on CIDA 
and the aerogel collector. As a result of the axially symmetric model, radial 
features on the night side of the nucleus appear (see figure ^) . Unfortunately, 
we can not prove or disprove the existence of these features by the Stardust 
measurements, because Stardust is passing by the day side of P/Wild 2. The 
prove or disprove of the existence of such features remains a task for future 
ground based observations and space missions. 

On the basis of the Galileo and Ulysses measurements we predict that during 
the interplanetary cruise phase 25 mainly small (m « 10"^^ g) particles will 
hit the CIDA target. As the Ulysses data indicate, Stardust's timing with 
the solar magnetic cycle is unfortunate (see section 3.3). To confirm that the 
decrease in the total fiux of interstellar particles measured by Ulysses is due to 
their interaction with the solar wind magnetic field, modeling of this interaction 
is under way [Q. Even if the total number of interstellar particles detected 
by CIDA could be low, the scientific significance of the data, in any case, is 
high because the composition of LIC dust is to be directly measured for the 
first time by CIDA. The reliability of this measurement depends very much 
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on the ability to identify interstellar impacts clearly. Interplanetary particles 
are a possible source of "contamination" for CIDA's measurements. Table |^ 
shows that 5 interplanetary dust particles with masses greater than 10~^^ g are 
expected to hit the CIDA target during the cruise phase. For the determination 
of this number we have used the five-population-model of interplanetary dust 
by Staubach et al. (1997) Q. So, interstellar particles should dominate the 
impacts on the CIDA target during cruise, but one has to make sure that an 
individual detection is really from an interstellar particle. 

We expect the aerogel collector to contain about 120 interstellar particles 
inside its interstellar collection layer. Of these, 40 should be large particles 
(size > 1 /im) that are probably the only ones which are extractable from the 
aerogel. Due to their high impact velocities, these large particles might be able 
to penetrate very deep into the aerogel and are possibly altered by the entry 
process. The distribution of impact angles can be determined by measuring 
the geometry of the entry-tracks. This will give valuable information on the 
dynamics of interstellar particles in the solar system. The highest amount of 
information is of course contained in the particles themselves when put into 
the laboratory and analyzed for their chemical and mineralogical properties. 
Like CIDA, the collector might also contain a contamination by interplanetary 
particles on both sides, the interstellar and the cometary dust collection layer. 
As shown in table || we predict 25 interplanetary particles with masses greater 
than 10"^'^ g to hit the cometary side of the collector and 20 to hit the interstellar 
side. 

If the interplanetary particles can be identified as such, it might be interest- 
ing to compare them with the IDPs collected by stratospheric aircraft M, p3] . 
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A Spatial Density in the Coma 

The number of particles which leave the coma per unit time and solid angle is 
given by Ze^i('&), where 1!) denotes the angle between ejection and Sun direc- 
tion. The velocity of the particles is given by Ve,i{^). For the computation we 
introduce a reference frame with the comet nucleus in its origin: The z-axis 
is pointing to the anti-solar direction. The observers position is inside the z- 
X plane. The position of a particle which was released in the (p) direction 
before time t with respect to this reference frame is given by (see figure : 

(cos{(p) sin(?9) \ 
- sm{tp) sin(z?) (26) 
-cos{d) J 

with r = tVe^^{'&) 

The phase angle of the observation is denoted by a. The projected position of 
a dust particle is given by (compare figure EOt): 




= P 



cos(?7) 
sin(77) , 

^ \ (27) 

- cos(q;) cos{(p) sin(i9) — sin(Q;) cos(i3) 

— sm{(p) sin('i?) 

The X-axis in the image plane is the projected anti-solar direction. Alternatively 
to the cartesian coordinates in the projected plane {x',y') a point in the image 
plane can be represented by the cylindrical coordinates (p, rf) , where the 77 = 
direction points towards the (projected) anti-solar direction. The number of 
particles which are ejected in the time interval dt and with solid angles inside 
dip sa\{'d)d'd is given by: 

dN = Ze4^)dtdif sm{i})di} 



drdipsmid)dd (28) 
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Next, we determine the contribution of a cone of given angle to the projected 
number density is computed. The total number density is then obtained by 
integrating over all cones in a second step. 

For computing the number of particles inside an infinitesimal area element 
in the projected plane dA = p dpdrj due to particles of the cone d, one needs to 
know for which radial distance r and azimuthal angle Lp there are contributions 
to the chosen point in the image plane {p,ri). This can be found by solving 
equation (|2^) for r and cos((p) (with constant -d). This leads to 

ri,2 = 27Qr~ 2rT (~ (^os{t]) sin(a) cos{d) 

cos^(v) — cos^(a) 

±cos{a)y/u^ -cos2(i?)^ (29) 
If sin(2a) cos(i9) . 2r \ 

± cos(r/) -cos2(z9)) , (30) 



where u = u[a,'q) = \J\ — (sin(77) sin(Q;))2. For u < |cos(i?)| the cone i} does 
not contribute to the point with angle 77 in the image plane. Furthermore, since 
r = Ve,i{'&)t > 0, only the solutions with r > have to be considered. The 
number of solutions for a given angle r], phase angle a, and cone •& is denoted 
by N{a, rj, ■&) and is given for different conditions in table ||. 

Finally, we calculate which infinitesimal azimuthal-radial element dr dip cor- 
responds to the chosen surface element dA = p dpdrj in the image plane. This 
relation is given by the determinant of the Jacobian matrix: 

dN = • det ( g ^\dpdTj sm{^)d§ 

\ dp dr) J ^21) 



Computation of the partial derivatives yields 

Since the last equation gives only the contribution to the number density due 
to a single solution of equations (^^ and (^^ the contribution of all solutions 
has to be summed up to find the total number of particles in dA = p dpdrj due 
to cone d: 

^=N{a,rj,i))^^^=2==sinm^ (33) 
dA pVe,i(v) ^u^ - cos2(i?) 

Finally, the projected number density is found by integrating over all cones 
which contribute to the chosen point (p, 77) : 

— (p,77) = / N{a,rj,-d)- 
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•^__=dcosW (34) 

For the calculation of the intensity which is received from the dust coma the 
number of particles inside a circular field of view with radius p is needed and 
can be found by integration of the projected number density over the field of 
view (note that (A^(a, r?, i?) + A^(a, tt - ??, z?))/2 = 1): 



'■2- Ze,,(7?) 



Jo J-u 



rfcos(^) dr] (35) 



= y' 




z 



Figure 20: We define the reference frame with the comet in its origin such that 
the z-axis points in the anti-solar direction and the observers position is located 
in the plane spanned by the x- and ^;-axis. In the projection plane, the y'-axis 
is parallel to the y axis and the x'-axis points in the direction of the projected 
anti-solar direction. Alternatively to the cartesian coordinates, a point in the 
projection plane can be represented by the cylindrical coordinates (p, 77), where 
p is the distance of the point from the projected comet position and r] the angle 
from the anti-solar direction. 



B Spatial Density of Hyperbolic Particles 

We derive the number density of particles, which initially enter the Solar System 
on parallel trajectories with velocity T4o and number density rioc, assuming 
that the particles are moving on Keplerian trajectories. For this purpose an 
infinitesimal volume element dV at distance r from the Sun and angle -d with 
respect to the direction the particles enter the Solar System is considered (see 
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Figure 21: The sketch shows the volume element which contains the particles 
originating from the volume element cIVq. The particles on the right boundary 
of the volume elements reach the left boundary in the time dtg ■ 



figure [21]) . The particles which cover the volume element dV has started in the 
time interval dta from outside the impact parameter interval db and covered at 
a large distance from the Sun the volume element dVo — 2'Kb dbVoodtQ. The 
number of particles inside both volume elements is given by dN = fioo dV^ . The 
size of the volume element dV can be computed by 

dv 

dV = 2-Kr sin(i5) dr Vi_ dto = 27rr —db V± dto. (36) 

Due to the preservation of the angular momentum on a Keplerian trajectory the 
velocity perpendicular to the radial direction can be expressed by V± ~ bVoc/r 
and the size of the volume element can be written 

dv 

dV = sm{^)—dVo. (37) 

Thus, only the derivative of the radial distance r with respect to the impact 
parameter needs to be computed. The radial distance r of a non radial Keplerian 
trajectory can be parameterized: 

P 



l + ecos(z?-?9o)' 



(38) 



where p, e and "do are constants which has to be expressed by the initial condi- 
tion: 



cos(??o) = 
sin(i?o) = 

d ^ ^, (39) 
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where /i = jMq{1 — j3) is the effective gravitational constant of the Sun and e 
is positive in the case /i > and negative for /i < 0. By inserting the constants 
p, e and do in equation (^8|), we find the radial distance as a function of the 
impact parameter h: 

9 1 

r(6) = r (40) 

^ ' 1 - cos{d) + I sin(i9) ^ ' 

By computing the derivative with h an expression for the number density at the 
position (r, d) is found: 

dN 5^ 

" = Tn? ^ ■ ^QMo^ woTT > (41) 

dV rsm(w)|2o — rsm(?y)| 

For computing the number density at a given position in the Solar System, the 
impact parameter 6 as a function of r and •& is needed, what can be found by 
solving equation (EO) for b: 



b = 



rsin(i9) //rsinft?) , 
— — ^ +rd(l-cos(^)) 



(42) 



For an attractive effective gravitational constant /i > 0, i.e. d > 0, every 
position can be reached by two trajectories (apart from the axes d — and 
i9 = 7r). For a repulsive effective gravitational constant /z < i.e. d < 
every position outside the parabola |rsin(i?)| — d^ + djl ■ rcos(i?) is reached 
by two different trajectories, no trajectories reach the interior of the parabola 
and every place on the parabola is reached by exactly one trajectory. Since 
equation ( p] ) gives only the contribution to the number density of one solution 
the contributions due to the different solutions which reach a point has to be 
added up to calculate the total number density at a point. On the axis i? = tt 
for yu > and on the parabola for /i < the volume element dV becomes zero. 
Therefore the number density becomes infinite at those positions. However, the 
infinities are of that kind that any integral of the number density over a finite 
volume is finite. 

Because the representation of a Keplerian orbit as it is given in equation 
( ^8| ) cannot be applied to radial trajectories i.e. impact parameter 6 = and 
I? = 0, the number density at the axis = cannot be computed directly by 
equation (^), but the formula must be smoothly continued to this axis: 



(43) 



The upper sign applies to the contribution to the number density due to the 
particles approaching the Sun and the lower sign to the departing particles. 
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